Ferrofluidity in a two-component dipolar Bose-Einstein Condensate 
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It is shown that the interface in a two- component Bose-Einstein condensate (BEC) in which 
one component exhibits a dipole-dipole interaction spontaneously forms patterns similar to those 
formed in a magnetic liquid subject to a magnetic field. A hexagonal pattern, hysteretic behavior, 
and soliton-like structure are numerically demonstrated. A phenomenon similar to the labyrinthine 
instability is also found. These phenomena may be realized using a ^^Cr BEC. The periodic density 
modulation in the superfluid ground state offers evidence of supersolidity. 

PACS numbers: 03.75.Mn, 03.75.Hh, 47.65.Cb 



When a magnetic liquid (a colloidal suspension of 
magnetic fine particles) is subjected to a magnetic field 
perpendicular to the surface, the liquid is magnetized 
and the surface undergoes spontaneous deformation into 
characteristic patterns shaped like 'horns' growing from 
the liquid. This surface instability is known as the 
normal-field or Rosensweig instability [1] and is the sub- 
ject of active research |2|, |3|, y, l5|, Isl, |7|, l8( . The system also 
exhibits a variety of phenomena, such as hysteretic be- 
haviors [1] , transition between hexagonal and square pat- 
terns [sf, and stabilization of a soliton-like structure [8|. 
Because of the visual appeal of the pattern formation, the 
dynamics of the magnetic-liquid surfaces are introduced 
even in art 

In the present Letter, we show that a Bose-Einstein 
condensate (BEC) of an atomic gas with a strong dipole- 
dipole interaction [13] exhibits instabilities and pattern 
formations similar to those in magnetic liquids. The sys- 
tem considered here is schematically illustrated in Fig. [H 
A two-component BEC is used, in which the atoms in 
component 1 have a magnetic dipole moment and the 
atoms in component 2 are nonmagnetic. A magnetic field 
is applied to fix the direction of the fully polarized mag- 
netic dipole (say, in the z direction). A magnetic-field 
gradient in the z direction pulls the atoms in component 
1 in the —z direction, and then the two components are 
phase-separated for a strong field gradient. The surface 
deformation in a magnetic liquid by the Rosensweig in- 
stability occurs when it lowers the sum of gravitational, 
surface, and magnetic energies [3]. In analogy with a 
magnetic-liquid system, the field gradient plays the role 
of gravity, and the quantum pressure and contact in- 
teratomic interactions create an interface energy, corre- 
sponding to the surface tension of a magnetic liquid. 

A two-component BEC is used because in order for 
the density pattern to be formed in a single-component 
BEC, the dipole-dipole interaction must overcome the 
contact interaction, which leads to dipolar collapse 
In fact, the parameter region in which a biconcave den- 
sity pattern is formed is in the immediate vicinity of the 
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FIG. 1: (Color) Schematic illustration of the system. A 
two-component BEC is confined in a trapping potential, in 
which component 1 has a magnetic dipole moment and com- 
ponent 2 does not. The magnetic dipole is polarized in the 
z direction by the magnetic field. The two components are 
phase-separated due to the field gradient dB/dz < 0. 



region of the collapse instability lH, T^. In Refs. 3, 1 
quasi-two dimensional (2D) systems are considered in or- 
der to avoid the collapse, and 2D solitons are shown to 
be stabilized by the dipolar interaction. Using a multi- 
component BEC is another possible way of realizing pat- 
tern formation without suffering a dipolar collapse, since 
texture formation is much easier than density-pattern 
formation. For example, spin domains can be formed 
by a small spin-dependent interaction [l^. 

In this Letter we show that a hexagonal pattern of 
density peaks emerges on the interface between nonmag- 
netic and fully polarized ^^Cr BECs [10]. We find that 
the number of peaks can be controlled by the strength of 
the field gradient. Hysteresis also occurs between flat and 
single-spike interfaces, which resembles soliton-like struc- 
tures in magnetic liquids [8]. We show that various pat- 
terns are formed, including a labyrinthine pattern [17j]. 

We consider a two-component BEC described by the 
macroscopic wave functions tpi and in the zero- 
temperature mean-field approximation. The wave func- 
tions obey the nonlocal Gross-Pitaevskii (GP) equations 
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given by 

+ J U{r-r')\^i{r')fdr'^^i, (la) 

= (-^ + V^2 + 522|^2p + 5i2|^ip)^2,(lb) 

where and Vn are atomic mass and trap potential 
for component n = 1 or 2. The interaction coefficients 
are given by gnn' = 27rfi /Mnn' with and Mnn' 
being the 5-wave scattering length and the reduced mass 
between components n and n^ The dipole-dipole interac- 
tion has the form [/(r) = /ioM^(l — 3z^/r^)/(47rr^), where 
/io is the permeability of vacuum and /i is the magnetic 
dipole moment of the atoms in component 1. For simplic- 
ity, we assume that the two components have the same 
mass M and the same number of atoms N/2 and that 
they experience the same axisymmetric harmonic poten- 
tial given by Vn = M/2[uj'^{x'^ + y'^) + ^'^z'^], where uj± 
and 00 z are radial and axial trap frequencies. Gravity 
only shifts the origin and can be neglected. The mag- 
netic field depends only on and dB/dz is uniform. We 
employ the mj = —3 magnetic sublevel of a ^^Cr atom 
for component 1, for which /i = 6/iB with fiB being the 
Bohr magneton and an :^ lOOae [18] with as being the 
Bohr radius. We assume an a22 — <^12- 

We first study stable states of the system, obtained 
by replacing i with —1 on the left-hand sides of Eq. ([T]). 
The numerical propagation is performed using the Crank- 
Nicolson scheme and the dipolar part is calculated us- 
ing a fast Fourier transform. The initial state of the 
imaginary-time propagation is a flat-interface state 
for a sufficiently large \dB/dz\ with small initial noise 
i/jn = + ^n, whcrc represents a small complex num- 
ber randomly chosen on each mesh. The small noise is 
added to break the axisymmetry of the system. 

Figure [2] (a) shows an isodensity surface of component 
1 for an = a22 = <^i2 = lOOaB and dB/dz = —600 
mG/cm. The hexagonal mountain-like pattern on the 
interface between the two components is similar to the 
Rosensweig pattern on a magnetic-liquid surface. Since 
the density is high and hence the dipole-dipole interac- 
tion is large at the center of the trap, the peaks around 
the center are higher than those on the periphery. Fig- 
ure [2] (b) shows the column-density profiles integrated 
along the z axis. It is interesting that the left and mid- 
dle panels of Fig. [2] (b) exhibit phase separation in the 
x-y plane, even though the scattering lengths do not sat- 
isfy the immiscible condition. This is because the density 
of component 1 is large at the peaks due to the dipole- 
dipole interaction, which repels component 2. In fact, 
the total density [right panel of Fig. [2] (b)] also has a 
hexagonal pattern. 

When the field gradient becomes tight, the number and 
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FIG. 2: (Color) Stable state of the two- component dipolar 
EEC for ail = ^22 = ai2 = lOOaB and /x = 6/jb in an axisym- 
metric potential with {uj±,uJz) = 27r x (100,800) Hz. The 
number of atoms is A/" = 4 x 10^ with an equal population 
in each component and the field gradient is dB/dz = —600 
mG/cm. (a) Isodensity surface of component 1. The color 
represents the z coordinate normalized by (?i/Muj±_Y^'^ . (b) 
Column densities f \ipi\'^dz, f \ip2\'^dz, and jQipil"^ -\-\ip2\'^)dz 
normalized by NM(jj±/h. The field of view is 50 x 50 /xm. 

height of the peaks decrease and eventually the interface 
becomes flat. The appearance and disappearance of the 
pattern exhibits hysteresis with respect to a change in the 
field gradient. Figure[3](a) depicts stable states obtained 
from the imaginary-time propagation of the GP equation 
as \dB/dz\ is increased or decreased. When \dB/dz\ is 
increased from a small value [red arrow in Fig. [3] (a)], the 
hexagonal pattern of the peaks changes into a single peak, 
which then disappears at dB/dz :^ —1054 mG/cm. On 
the other hand, when \dB/dz\ is decreased from a large 
value [blue arrow in Fig. [3] (a)], the interface remains 
flat until dB/dz —1028 mG/cm. Therefore, there is a 
bistable region at \dB/dz\ 1028-1054 mG/cm, in which 
the single-peak state and the flat-interface state are both 
stable. In the bistability region, the energies of the two 
states are almost degenerate. Such hysteretic behavior 
has been predicted [3^, J} and observed [4, S] for magnetic 
liquids, in which the applied magnetic field and hence the 
magnetization is changed. For a wider pancake-shaped 
BEG, bistability between the flat-surface and hexagonal- 
peak states may be observed as in a magnetic liquid. 

The single-peak structure in the bistable region in 
Fig. [3] (a) is similar to the "ferrosoliton" [8] in a magnetic 
liquid, which shows a stable soliton-like peak. Figure [3] 
(b) shows the density profile of the cross section of the 
single-peak state at dB/dz :^ —1040 mG/cm. Around 
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FIG. 3: (Color) (a) Hysteresis with respect to a change in the 
field gradient. The red and blue arrows indicate that a single- 
peak or flat interface is stable for \dB/dz\ ^ 1054 mG/cm 
or \dB/dz\ ^ 1028 mG/cm. Beyond these thresholds, the 
two stable states convert to each other (dotted arrows). The 
isodensity surfaces are shown for dB/dz :^ —900 and —1040 
mG/cm. (b) Density profile of the cross section on the y — ^ 
plane for the single-peak state at dB/dz ^ —1040 mG/cm. 
The field of view is 44 x 6.7 fxm and the density is normalized 
by N{Mou±/?i)^^'^ . The parameters are the same as those in 
Fig. [2] except for dB/dz. 
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FIG. 4: (Color) Isodensity surface (upper figures) and col- 
umn density (lower panels) of component 1 for various stable 
states. The difference in (a)-(c) is only in the initial perturba- 
tions in the imaginary-time propagation. The field gradient is 
dB/dz — —200 mG/cm. The other parameters and the color 
scales are the same as those in Fig. [2] 



the peak, the interface oscillates in a concentric manner, 
as in the ferrosoliton in a magnetic liquid 

To investigate various met ast able patterns, we apply 
initial perturbations to trigger the pattern formation. 
Figure [4] (a) is obtained by applying an additional mag- 
netic field (X sm.(kxX^kyy) for t < 1 ms in the imaginary- 
time propagation, where kxd^ = and kya^ = 1 with 
a±_ = {%/Mu^Y^'^ . We find that the imaginary-time 



FIG. 5: (Color) (a) Isodensity surface of component 1 and 
column density of (b) component 1 and (c) component 2 for 
ail — lOOttB, ai2 — a22 — 120aB, dB/dz — 0, and N — 
8 X 10^. The other conditions and the color scales are the 
same as those in Fig. [2l 



propagation relaxes the system to the stripe pattern as 
shown in Fig. [4] (a). This pattern is robust against small 
additional noise. The initial temporary magnetic field for 
Fig. [4] (b) is (X cos(y^x^ + y'^ /ai.)- The concentric pat- 
tern in Fig. [4] (b) is stable against axisymmetry-breaking 
noise. Figure [4] (c) is obtained without an additional 
magnetic field, with only a small noise added to the 
initial state. These results imply that there are many 
metastable states with various patterns. In Fig. [H (c), 
the peaks partly merge into each other. The merging of 
the peaks tends to occur for a large dipole interaction 
and small field gradient. For dB/dz = —600 mG/cm, 
the stripe and concentric patterns in Figs. [4] (a) and [H 
(b) are unstable against forming peaks. Although the 
square lattice of the peaks is found to be unstable for the 
parameters in Figs. [2] and [H a hexagonal-square-lattice 
transition may occur for other parameters, as in magnetic 
liquids 0,[5|. 

Figure [5] shows a stable state for dB/dz = which has 
an intricate pattern similar to the labyrinthine pattern in 
a thin layer of magnetic liquid confined with an immis- 
cible nonmagnetic liquid [l7|. The relation between the 
labyrinthine instability in a magnetic liquid [2] and that 
in the present system remains to be clarified. 

Next we study the dynamics of the pattern formation 
in a nondissipative system by solving the real-time prop- 
agation of Eq. ([T]) . Figure [6] shows the time evolution of 
the system, where the initial state is a flat-interface state 
for dB/dz = —1200 mG/cm, as shown in Fig. [6] (a). In 
order to break the axisymmetry, a small initial noise is 
added to the initial state. We let the system evolve in 
time with dB/dz = —600 mG/cm, for which the stable 
state has hexagonal peaks, as shown in Fig. [2] (a). The 
density fluctuations grow into many peaks, as in Figs. [2] 
(b) and[2](c), due to the Rosensweig instability. We note 
that the characteristic wavelength in the pattern at the 
early stage [Fig. [2] (b)] is smaller than that in the later 
pattern [Fig. [2] (c)]. This difference indicates that the 
most unstable wavelength in the linear regime does not 
correspond to the final stable pattern determined by the 
nonlinear interaction. A similar situation also occurs in 
magnetic liquids [6] . If the phenomenological dissipation 
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(a) t = ms (b) t = 5 ms (c) t = 40 ms 




FIG. 6: (Color) Time evolution of the isodensity surface 
(upper figures) and column density (lower panels) of com- 
ponent 1. The initial state is the flat- interface state for 
dB/dz = —1200 mG/cm and the field gradient changes to 
dB/dz = -600 mG/cm at t = 0. The field of view is 50 x 50 
fim. The parameters and color scales are the same as those 
in Fig. El 

is taken into account by replacing i with i — 0.03 [19] in 
Eq. ([T]), the hexagonal pattern is formed at ~ 50 ms. 

A candidate for component 2 is the ^83 mj = state 
of a ^^Cr atom. Though the mj 7^ —3 states are unstable 
against dipolar collisions, the lifetime is long enough (~ 
s) [10] to observe the pattern-formation dynamics. The 
spin exchange dynamics, such as mj = + ^ —3 + 3, 
can be suppressed using, e.g., the microwave-induced 
quadratic Zeeman effect [21]. The scattering lengths 
measured in Ref. [18i] give a22 — 60.7aB + ao/7 and 
ai2 — 67.8aB, with ao being the scattering length for the 
colliding channel with total spin 0, which is unknown. 
The scattering length an can be controlled using the 
magnetic Feshbach resonance [2o|. The inhomogeneity 
of the magnetic field ~ IG/cm x 10 /im is much smaller 
than the resonance width 1.7 G [13] • We can therefore 
obtain an ~ a22 ~ <^12 ~ OOae by changing an, if |ao| is 
not very large. We have confirmed that the Rosensweig 
pattern emerges for these scattering lengths. Another 
possibility for component 2 is an alkali atom, for which, 
however, ai2 has not been measured yet. The pattern 
formation for various scattering lengths merits further 
study. 

An important difference between the present system 
and magnetic liquids is that the present system is a su- 
perfluid, while a conventional magnetic liquid is a normal 
fluid. In the Rosensweig pattern, the total density shows 
a periodic modulation (i.e., diagonal order), while the 
present system by construction possesses an off-diagonal 
order, presenting strong evidence of super solidity. An- 
other difference is that the present system is fully po- 
larized irrespective of the strength of the magnetic field, 
while the magnetic liquid is paramagnetic. The difference 
between a compressible gas and an incompressible liquid 
is also important. It is remarkable that both systems 
exhibit similar phenomena despite these differences. 

In conclusion, we have studied pattern formation on 



the interface in a two-component BEC, in which one com- 
ponent exhibits a dipole-dipole interaction. We found a 
rich variety of inter facial patterns, including Rosensweig 
hexagonal peaks and a labyrinthine pattern. We also 
observed hysteretic behavior and the ferrosoliton, as in 
magnetic liquids. We expect that the two-component sys- 
tem proposed here will provide a new insight into surface 
and interface physics with long-range interactions. 
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